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Abstract 

The systems exhibiting quantum phase transitions (QPT) are investigated 
within the Ising model in the transverse field and Heisenberg model with easy- 
plane single-site anisotropy. Near QPT a correspondence between parameters 
of these models and of quantum (ft 4 model is established. A scaling analysis 
is performed for the ground-state properties. The influence of the external 
longitudinal magnetic field on the ground-state properties is investigated, and 
the corresponding magnetic susceptibility is calculated. Finite-temperature 
properties are considered with the use of the scaling analysis for the effective 
classical model proposed by Sachdev. Analytical results for the ordering tem- 
perature and temperature dependences of the magnetization and energy gap 
are obtained in the case of a small ground-state moment. The forms of depen- 
dences of observable quantities on the bare splitting (or magnetic field) and 
renormalized splitting turn out to be different. A comparison with numerical 
calculations and experimental data on systems demonstrating magnetic and 
structural transitions (e.g., into singlet state) is performed. 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



The interest in quantum models of anisotropic spin and pseudospin systems is connected 
with that they describe miscellaneous magnetic and structural transitions. Examples of such 
transitions are transitions into singlet magnetic state in TbSb, Pr, P^Tl (see Ref. Q and 
references therein), NiSi 2 F 6 (see, e.g., Ref. 0]), and orientational and metamagnetic phase 
transitions under magnetic field 

The simplest model for the systems demonstrating a ground-state quantum phase transi- 
tion (QPT) is the Ising model in the transverse field. This model is convenient for description 
of structural transitions in quantum crystals [§-0]. It can be also applied to describe mag- 
netic systems where both lowest and next energy levels are singlets. A more complicated 
first-principle model for spin systems in a strong crystal field is the Heisenberg model with 
an easy-plane single-site anisotropy; it is applicable in the case where next-to-lowest energy 
level is doublet. As well as transverse-field Ising model, this model also demonstrates QPT 
(the ground-state magnetization vanishes with increasing the anisotropy parameter). 

A number of approximate methods were applied to study the transverse-field Ising model 
P |16| at d > 1 and the Heisenberg model with easy-axis anisotropy [17— 21 1 . However, all 
these methods (except for numerical ones) are applicable only not too close to QPT. In 
particular, they lead to Gaussian values of the QPT critical exponents. Thus analytical 
consideration of the ground-state and finite-temperature properties in the vicinity of QPT 
is still an open problem. 

The case, where the system is close to QPT, is characterized by a small ground-state 
moment and low transition temperature. Such a situation is reminiscent of weak itinerant 



magnets [22 1 . An analysis of the ground-state QPT was performed in Refs. [|3|,[2J] . It was 
shown that the upper critical dimensionality for such transitions is d^ = 4 — z with z being 
the dynamical critical exponent. This conclusion has a general character. In the present 
paper we consider only systems with z = 1, which holds for the transverse-field Ising model 
and anisotropy-induced QPT in the Heisenberg model. (Note that this is not the case for 
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the QPT induced by magnetic field in degenerate systems with n > 1 component order 



parameter since here z — 2, see Ref. [0). Thus for d > 3 the critical exponents are the 
Gaussian ones, while for d < 3 they deviate from the corresponding mean-field values and 
can be calculated with the use of the 3 — e expansion. For the critical dimensionality d = 3 
the ground state properties contain logarithmic corrections. 

Sachdev |26j proposed a three-stage method of treating finite-temperature properties of 
the systems near QPT. At the first stage, ground-state renormalizations are performed. At 
the second stage, the non-zero Matsubara frequencies are integrated out to obtain an effective 
classical action. Finally, perturbation theory for the effective classical model is applied. This 
method ensures correct analytical properties of the resulting theory. While ground-state 
renormalizations are non- universal, finite-temperature properties, being expressed through 
quantum-renormalized ground-state parameters, turn out to be universal. 

The approach of Ref. 0] is based on a continuum model, namely, the quantum 4 
model. This model is sufficient to express finite-temperature properties near QPT through 
the non-universal ground-state properties, but insufficient to obtain correct results for the 
latter. A convenient method to consider the lattice spin systems near their critical dimen- 
sionality is the expansion in the formal quasiclassical parameter 1/S. Its applicability is 
connected with the fact that near d c the effective interaction of spin waves is small (except 
for a narrow critical region where the e-expansion can be easily developed to correct the 
description of the critical behavior). For the Heisenberg model such a situation occurs for 
temperature transition near the lower critical dimensionality d~ = 2. This provides suc- 
cess of the renormalization-group (RG) approach for the description of thermodynamics of 
d = 2 (Ref. |p7| ) and d = 2 + e (Refs. [p8p7| ) Heisenberg magnets, and also quasi-2D and 



anisotropic 2D magnets |29| not too close to T c . 

For QPT in highly-anisotropic spin systems the l/S'-expansion works well near the upper 
critical dimensionality d£ = 3. In this case there are excitations, which are almost gapless 
near QPT (they are analogous to the spin- wave excitations in Heisenberg magnets). Besides 
that, for the ordered degenerate systems (with n > 2) there are always Goldstone modes 



with zero energy gap and the l/S'-expansion becomes applicable at arbitrary anisotropy 
below its critical value. The situation is more complicated for finite temperatures, since 
close to temperature transition the system behaves as a corresponding classical magnet and 
therefore the picture of excitation spectrum differs from that at T = 0. 

The aim of the present paper is to apply the above- discussed concepts for calculat- 
ing ground-state and finite-temperature properties of transverse- field Ising model (n — 1) 
and Heisenberg model with strong easy-plane anisotropy (n = 2). To this end we apply 
perturbation theory (which is in fact an expansion in 1/S) to the original lattice models 
(not to their continuum analogs), which enables us to calculate non-universal ground-state 
quantum renormalizations. After that we combine perturbation results for short-wave fluc- 
tuations with the results of the 3 — e RG approach for the long-wave fluctuations to correct 
the results of perturbation theory. Finally, we consider finite-temperature properties within 
RG approach for the effective continuum classical model. 

The plan of the paper is as follows. In Sect. 2 we discuss the Ising model in the transverse 
field. We consider corresponding mean-field results, construct the perturbation theory in 
1/S and apply a scaling approach to investigate ground-state and thermodynamic properties, 
in particular the influence of external magnetic field. In Sect. 3 the Heisenberg model with 
easy-plane anisotropy is considered in a similar way. In Sect. 4 we discuss the results obtained 
and compare them with experimental data on systems exhibiting structural and magnetic 
transitions. Some details of calculations are presented in Appendices. 

II. TRANSVERSE-FIELD ISING MODEL 
A. The formulation of the model and the mean-field approximation 

We consider the Hamiltonian of the Ising model in the transverse field 

w = ~Es?s?- n E<s" (i) 
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where / is the exchange parameter. This model can describe singlet magnetic systems. 
A derivation of such a model for Heisenberg magnets with strong single-site anisotropy is 
presented in Appendix A. The model (|I|) describes also structural transition in quantum 
crystals (cooperative Jahn- Teller effect, see Ref. J7|) where two lowest energy levels are 
singlets. In this case I = A 2 , = A 1; where A 12 is the energy-level splitting at T = 
and T > T c respectively. For further purposes it will be useful to consider the model ([!]) for 
arbitrary values of (pseudo-) spin S. 

At Q = the model (H) coincides with the Ising model and thus the order parameter 
S = (S x ) — S in the ground state. With increasing Q, the model ([I]) demonstrates a 
quantum phase transition where S vanishes. The one-dimensional S = 1/2 transverse-field 
Ising model in the ground state can be solved rigorously [[3(J . In particular, it can be reduced 
to the two-dimensional Ising problem at finite temperatures fl31"l , so that critical exponents 



for both the phase transitions coincide. The transverse-field Ising model with d > 1 requires 
approximate methods. 

The mean-field theory |6]|| yields the critical field Qq = IqS, and the equation for the 
order parameter at Q < Qq reads 

Q 



B s (HJT) = 1 (2) 



where 



B s (x) = (1 + 1/2S) coth(l + l/2S)x - (1/2S) coth(x/2S), 
B 1/2 (x) = (1/2) tanh(x/2) (3) 

is the Brillouin function, H e = (fi 2 + tt 2 Q S 2 / S 2 ) 1/2 , J = 2Id. Owing to the field fi, the value 
of (S z ) is finite in both ordered and disordered phase and reads 

(S z ) = -j^B s (HJT) . (4) 

It should be noted that at Q < Qq we have simply (S z ) = SQ/Qq- The critical temperature 
where S vanishes is determined for the physically important case S = 1/2 by 



2tanh- 1 (fi/fi ) ln[2/(l - fi/n )] 
(the last approximation is valid for 1 — Q/Qq 1)- Thus the mean- field theory predicts a 
very weak inverse-logarithmic dependence for the critical temperature near QPT in arbitrary 
dimensionality. This contradicts to the results of the scaling approach p^ , |26| both above 
and below the upper critical dimensionality cf+ = 3. 

To improve the mean-field theory, one has to take into account the collective excitations 
which are analogous to spin-wave excitations in Heisenberg magnets. The spectrum of these 
excitations in the random-phase approximation has the form |J 

E^nin-i^ + iZs 2 (6) 

in both ordered and disordered phases. Near QPT (at S(T = 0) C 1), these excitations 
become almost gapless and give dominant contributions to physical properties. 

The result of account of the collective excitations to first order in 1/JZ (where JZ is the 
radius of exchange interaction) |Tj| for d = 3 reads 



T c ~n 3/ yi-n/n . (7) 

This has a correct square- root behavior (see, e.g., [PP]). However, the logarithmic corrections, 
that occur for d = 3, are not reproduced by the result (|?p. Besides that, the 1/JZ expansion 
does not enable one to determine correctly the coefficient in (|7|) for not too large JZ. 

Another approach used in O] is to consider the excitations (|6|) self-consistently within 
the random-phase approximation (RPA) decoupling scheme for the sequence of equations of 
motion. Unlike the 1/JZ expansion, this procedure gives a possibility to take into account 
the reaction of the RPA excitation spectrum ([]) on the deviation of the critical field from Q . 
However, corrections to mean-field ground-state parameters turn out to be small enough, 
and at finite temperatures the RPA magnetization shows a double- value behavior with first- 
order temperature phase transition. Authors of Ref. JO]] consider also a generalization of 
RPA, the two-site self-consistent approximation (TSCA) which gives a possibility to include 



partially correlation effects. This approximation gives more satisfactory results than RPA. 
However, it predicts first-order character not only for the temperature transition, but also 
for QPT. 

One should mention also the papers |p~2|JT3[1 where high-temperature series expansions 
(HTSE) and ground-state perturbation theory (GSPT) were used. Although these expan- 
sions gives consistent results for the critical field, their applicability near QPT is ques- 
tionable. Recently some results of GSPT and HTSE have been confirmed by numerical 
correlated-basis-function analysis [[l(| . 



Below we use the 1/5* expansion to treat ground-state and finite temperature properties 
of the transverse- field Ising model. Unlike the 1/1Z expansion, it takes into account the 
"spin-wave" excitations already in zeroth order of perturbation theory. Contrary to RPA 



TT| , this is a systematic expansion, and therefore scaling corrections can be easily calculated. 
It should be also noted that the 1/S expansion differs from the ground-state perturbation 
theory used in Ref. [|13j where the expansions in powers of Q/I and I/Q are used for the 
ordered and disordered phases respectively. Indeed, the 1/S" expansion treats both the terms 
in the Hamiltonian ([I]) on equal footing and thus yields physically correct results already in 
the first order in 1/S. 

B. Ground-state properties within the 1/S perturbation theory 

To construct perturbation expansion in a convenient form we use the spin coherent state 
approach [^|. The partition function is presented in terms of a path integral, 



Z = J D7rexp[-(S dyn +S st )}, (8) 



where 



l/T 

Sdyn — i£>y^ j J dr (1 — costV 





l/T 

S st = - J dr 
o 



dr 

IS 2 



2 (ij) 



(9) 



are static and dynamic parts of the action, = {n xi , TT yi , n zi } is three-component vector field 
with n 2 = 1 + 1 /S,&i and (fi are the polar and azimuthal angles of 7Tj in an arbitrarily chosen 
coordinate system (which does not need to coincide with the 7i x -7i y -n z coordinate system). 
Further we additionally rotate the coordinate system through the angle 9 determined by 

sin0 = (^>/(kl> (10) 

around ^-axis (in the disordered phase 9 = and the rotated coordinate system coincides 
with the original one). Then (tt x ) = in the new coordinate system. 

The calculation of two-point vertex function r(q, uj) of the fields n x , n y (the tilde sign 
is referred to rotated coordinate system), which is connected with matrix Green's function 
G of these fields by the relation r(q, uj) = G _1 (q, uj), is performed for both ordered and 
disordered phases in Appendix B and yields to first order in 1/S the result 



r±(q, uj n ) 



S 2 (I -I n + I A 2 ± ) iSuj n (w s + X /2 + Y /2)^ 



(11) 



^ iSuj n (w s + X /2 + Y /2) I S 2 D ± ) 

where w s = (l + l/2S)-\ 



q y 1 — Iq/ lo q 
A± and D± are the dimensionless temperature-dependent energy gap and the renormaliza- 
tion factor for the exchange parameter in the disordered and ordered phases respectively, 
their concrete expressions being specified below. The matrix static uniform spin suscepti- 
bility in the rotated coordinate system is expressed through Y as 

^ = ^(0,0) (13) 

where i,j = x,y. The renormalized spectrum of "spin-wave" excitations is determined by 
the condition detT(q, —iE^) = and to first order in 1/S has the form 



E n = S [1 + 1/2S - (X + Y )/2] yV^/o-Zq + ZoAl). (14) 
The quantum-renormalized critical field Q c is given by 



9l = \ + — 1 V 2io + / q 



Last two terms in this expression yield the first-order l/S'-correction to mean-field value of 
Q c . For 5=1/2 numerical calculation of integral in ([H^) yields the result Q c = 2. 44/ in 
the 3D case and fl c = 1.10/ in the 2D case. Thus the critical field is strongly renormalized 
by quantum fluctuations both in 3D and 2D cases. The critical field values obtained are 



considerably smaller than the corresponding RPA results [11], f2 c = 2.88/ and Vl c = 1.83/ 



and somewhat smaller than those obtained by HTSE [12| and GSPT |Tj|, fi c = 2.58/ 
and Q c = 1.54/. This demonstrates that considered first-order 1/S perturbation theory 
overestimates effects of quantum fluctuations (especially in the 2D case), but treats these 
fluctuations more correctly than RPA. 

In the disordered phase with (ir x ) = (Q > Q c ) the expressions for the ground-state 
energy gap and factor D + has the form 

1 — £+ 

A+{t)= JL E f 2/o + J q (l +t ) _ 27 + J q 1 

4S '«l/»-',(i-t)] •JWo - Q J 

and 

D + (t + ) = 1+ ^_~ X ° [1 + UA + (U)} (17) 



where 



and 



In the ordered phase (f2 < f2 c ) we obtain 

A 2 (t_,0) = t_(l-X^) [1 + A_(f_)], 
A_(t) = -2(1 -t)A + (t) 



t+ = i- n e /n. (is) 



l-t [2/ + / q (l + t)F (lg) 



85 V/ 1/2 [/o-(l-t)/ q ] 3 / 2 
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and 

D_ = 1 + Y - X (20) 

where 

t- = i - (o/q c ) 2 . (21) 

Consider now the observable quantities. The expression for the order parameter 
S(t-,T) = S(tt x ) at T = 0, Q < Q c reads 

S(t_, 0) = St 1 ! 2 [1 + 5(i_)] 1/2 [1 + 1/25 - (X + Y )/2] , (22) 
B{t) = -2(1 -t)A+(t) 

For the longitudinal susceptibility we have 

x ™ = cos 2 0x xx + sin 9 cos (jf* + X™) + sin 2 (23) 

where the tilde sign is referred to susceptibilities in the rotated coordinate system (recall 
that for the disordered phase 8 = 0). For the ordered phase the first summand in ( p3|) gives 
dominant contribution near QPT, and using the relation (|1^) yields in both the ordered 
phase near QPT and disordered phase the expression for the ground-state spin susceptibility 
through the gap in the excitation spectrum 

*" = TZmsry (24) 

In the limiting case of zero transverse field we have the trivial result S = S, and the 
excitation spectrum reduces to its mean-field form, E n = Q . At very large Q ^> Q c we 
reproduce again the mean-field result Eq = Q. This is a consequence of the fact that in both 
the limits Q = and Q — > oo quantum fluctuations are absent. Thus the 1 / 5-expansion gives 
a possibility to describe the ground-state properties for an arbitrary Q > 0. However, as we 
shall see below, some difficulties arise in the region Q w Q c where the quantum fluctuations 
are strong enough to modify considerably the results. A more detailed consideration of this 
region will be performed in Sect JII D| . 
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C. Influence of longitudinal magnetic field 



To consider the influence of the external magnetic field we add to the Hamiltonian the 
term 

AH = -HY,S?- (25) 

i 

The longitudinal magnetic field results in the appearance of nonzero (Sf) at any Q/I. The 
influence of both transverse and external longitudinal fields is, of course, equivalent to ap- 
plying one effective field which has the value (H 2 + VL 2 ) 1 / 2 and makes the angle arctan(i//fi) 
with the 7r x -axis. However, it is useful to consider these fields as two independent ones. 
Performing the calculations which are similar to those of Sect . [II B| and Appendix B we 



obtain to first order in 1/S the equation for the angle 9 of coordinate system rotation 

n-n r(6)cose - HcotO = 0, (26) 



where 



r(e) _ 1 , 1 1 2(I + / q Mfl)-I q cos 2 fl 

2S 45 q Jj ^)[/o^)-/ q OOB 8 e]' 



with (p(0) = sin 9 + (fl/Q c ) cos9. For a general Vt/VL C the solution of this equation is rather 
cumbersome. However, near QPT (i.e. at 1 — Q/Q c <C 1), where the angle 9 is small, one 
can expand (P5[) in 9 to obtain 

9 + H/2[I Sr(0)-Q] 9 H ^9 
9 = <( (28) 

9 H + 25r/39 H 9 < 6 H < 1 



where fl = V 2 ! 1 ~ ^M=), ^ = (2#/fi c ) 1/3 and 5r = r(^) - r(0). 
For the magnetization we derive 

_ = s(h = o) + ^ « i - n/flc, 

+ A - ^) ^ I 1 + > 1 " n M « °h < 1, 



where is determined by (P4|), and 
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B'K) 



E 



65^ q 



2(/ + / q )(l + ^/2) 



iq(l 



/ /o[/o(l + ^)-/ q (l-^/2)] 
The ground-state energy gap is given by 



2/p + I q 



(30) 



A 2 _(H = 0), 



0% < i-n/n c 



3/)2 



(31) 



^[1 + M)]> l-0/fi c <^<l, 



where 



(/o + 2/ q ) 



(32) 



12 5V/o 1/2 [/o(l + ^/2)-/ q (l-^)]3/ 2 ' 
The longitudinal susceptibility in the presence of magnetic field is still determined by 
(f23|), and again the first term gives main contribution near QPT. Alternatively, the same 
result can be obtained by direct differentiation of S [which is given by (|29"D1 with respect to 
H. 



D. Ground-state renormalizations near QPT 

The results of 1/S expansion can be applied only at not too small t±. Indeed, at d < 3 
the functions A and B contain terms which are divergent at t± — > as t± 3 ^ 2 (at d = 3, 
logarithmic divergences are present). The same situation takes place for the functions A' 
and B' which are divergent as O'jf 3 at 9h 0. Thus ane = 3 — d expansion can be developed 
within the RG approach to treat these divergences more correctly and to improve thereby 
the behavior of A± and S near QPT. Further consideration of this section is related to 
the critical region |1 — fi/fi c | <C 1. However, as it will be clear below, the results can be 
extrapolated to arbitrary Q, since in the limits Q <C Q c and Q ^> Q c they are smoothly 
joined with the results of the 1/S* expansion of section [11 B[ 

First we pick up the nonuniversal factors from S, A± by introducing the quantities 

S R (t, T) = [1 + 1/2S - (X + Y^/2}- 1 S(t, T)/S 
A ±R (t,T) = (1 - X^r 1 A±(t,T) (33) 
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Consider the continuum limit of the above theory. The action S = Sd yn + <5>st m this limit 
takes the form 

c/T 



Scant = ~ / <l' l r I (I T 



2m x (dn y /dT) + ?fj + (V^) 2 + m 2 7r 2 x 





c/T 

M r -id^, f j_~4 



+ -ld d rldrn x , (34) 



o 

,2 



where the parameters u, m and c, determined in such a way, are given by 



^cont 




Ccont 


— c 0, 


^cont 


= 2i+d, 



(35) 

c = (2d) 1 ^ 2 IS being the bare spin-wave velocity, 7? 2 = (I S 2 / Cq)ti 2 , tc 2 = (IS 2 /cq)ti 2 and the 
factor £ (£ = 1 — £_ in the ordered phase and £ = 1 in the disordered phase) is introduced 
to extend the region of applicability of results obtained to arbitrary t±. Note that the 
coefficients at first three terms of the quadratic part of (33) can be always chosen equal to 



their values in (|34D by appropriate rescaling of ir x , y and r. The model fl34D is completely 
equivalent to the quantum 4 model. Indeed, integrating out the field ir y we obtain 



1 



c/T c/T 



5 CO nt = - I d r I dr 



u 



(dn x y + m 2 n 2 x \+^J d d r J dr^ x (36) 
o " " o 



The continuum representation (EMI) determines the way in which the original lattice model 



can be renormalized. Following to the standard procedure |33[ we introduce the renormal- 
ization factors Zf for the ground-state parameters in the disordered and ordered phases by 
relations 

= Z x TT xR , TTy = Zy-KyR 

t± = (Zt/Z) t ±R , g = (Zf/Z 2 )g R (37) 
where the indices R denote quantum-renormalized quantities, 
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g = K 4 _ £ L e n e u cont (38) 

is the coupling constant, fi is a parameter with the dimensionality of inverse length, = 
[2 d - 1 7 r d / 2 T(d/2)}~\ and the factor L £ = T(l + e/2)T(l - e/2) [T(z) is the Euler gamma 
function] ensures the applicability of the one- loop order results for not small e [33]. For 
further treatment it is useful to represent the renormalization factors as 

zf = zUg)zr\g R ^) (39) 

where Z^ ont are the corresponding factors for the continuum model (|34]) that contain diver- 
gent terms (which are independent of lattice structure etc.) and Zu all the others (lattice 
dependent) corrections. It is important that the factors Z^ do not contain divergences. 



The expressions for Z-factors in the continuum model fl34]) are well known (see, e.g., Ref. 
). We use the cutoff scheme with cutting integrals over quasimomentum at A. Then to 
one-loop order we have 



2e V A 



zr = i + 3 -§{ 1 -^)- (4°) 

For our purposes it is convenient to put A = (2c?) 1 / 2 (the lattice constant is assumed to be 
equal to unity), rather than to pass to the limit A — ► oo (as it is usual in the quantum 
field theory). The expressions for Zu can be deduced by comparison the above results 
of perturbation theory for the original lattice model (Sections [II B| and |I C|) and standard 



perturbation results for the continuum model (|34]), see Ref. []33|]. We obtain 

= Z^y = 1 



r±X-l , , 9 ( 1 



[z&y 1 = i + At(t±) + 



2s \ t f 



{Zl,)- 1 = 1 + - + | (41) 
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Note that, unlike the factors Zf° nt , the quantities (f4T|) are defined only for integer e. As 
follows from fl39|), the determination of factors Zn enables one to consider the continuum 
model (|34]) with the parameters 



m 



cont 



c (1 + 1/2S - X 



(42) 



instead of the original lattice one. Thus the factors Zn represent the corrections owing to 
passing from the cutoff scheme in the original lattice model to that in the continuum model, 



cf. Ref. IE7I 



The flow functions for the coupling constant and energy gap have the standard form |33 

dg R 3 



i{9r) = v 



dp 
d\nZ™ nt 

dp 



£9R + ^9r 



~2 9r 



(43) 



The effective-Hamiltonian parameters g p , t p at the scale p! = pp as determined by these 
flow functions read 



9 P 



t, 



i -i 



P £ 9R 



-1/3 



(44) 



where g* = 2e/3 is the stable fixed point to one- loop order. 

We start the scaling procedure at p = A and stop it at p! = At l ± 2 (thus p = t± 2 ). For 
A± and S we obtain the results 

1 t± 



Aifl(t±,0) 



7 ± 



i+(3W2 £ )(i/4 /2 -i; 



1/3 



S R (t-,0) = t 1 J* 
where, according to ([37]), (fSJ), 



J LA 



U 



L2 



l + — -1 



1/3 



(45) 
(46) 



g R = (Zt L )-\2d)-^ 2 K^ e L e 



(47) 
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In the 3D case 



- 1 



ilni 
2 t 



(48) 



e Vt £ / 2 

so that the QPT critical exponents for the order parameter and the gap (inverse correlation 
length) are Gaussian one, 

(3 = 1/2, v = 1/2, (49) 

and logarithmic corrections are present. At the same time, in the 2D case we obtain 

(3 = 1/3, v = 7/12 (50) 

which are standard one-loop results for the one-component 4 theory in d+1 = 3 dimensions. 
In a strong enough longitudinal magnetic field (1 — Q/£l c <C 9 2 H <C 1) we obtain 

2 1 



Ai R (^,o) 



3Z'r 



L2 



Sr(9 h ,0) = 9 h < 



Z' 



l + (3W2e)(l/^-l)] V3 

1/3 



L4 



\2i 



L2 



1 + ^ ^-1 

2 £ ^ y 



(51) 



where 



J L2 



2e V 01 



(52) 



Thus, as well as for the dependences of ground-state properties on t, in the 3D case one has 
the mean-field value 



5 = 3 



(53) 



and the logarithmic corrections are present. At d = 2 we obtain the critical exponent 



5 = 9/2 



(54) 



Note that the scaling relations at d = 2 are slightly violated since the corresponding value of 
e is in fact not small and the e-expansion is applicable with a poor accuracy. However, this 
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violation is not too large (the value 5 = 5 can be calculated taking into account that the 
critical exponent 77 = to one- loop order), which indicates that the one-loop approximation 
gives adequate results even in this case. 

The ground-state parameters at zero longitudinal magnetic field are shown and compared 
with RPA and GSPT results in Fig.l for the 3D case and in Fig. 2 for the 2D case. The 
corrected values of Q c obtained from GSPT (see above), instead of those from first-order 1/S 
expansion, are used in the calculations. The l/S'-results for fl c are marked by arrows. One 
can see that, unlike results of RPA and l/S'-expansions, RG results have a correct critical 
behavior with critical exponents given by (^) and (0); besides that, they are very close to 
GSPT result for d = 3. For d = 2 the difference between RG' and GSPT results increases, 
which demonstrates that the ^-expansion has a poor accuracy here. Far from the quantum 
phase transition, the RG results coincide with those of 1/S perturbation theory. 

E. Finite-temperature properties near QPT 

At finite temperature the situation is more complicated, since not only "spin-wave" 
excitations, considered in previous sections, contribute to thermodynamic properties. At 
Q/Q c <C 1 we have T c ~ IS 2 and the phase transition occurs due to vanishing of ( | vr | } . 
The dominant excitations in this case are domain walls. Another situation occurs near QPT 
(1 — Q/Q c <C 1) where the temperature phase transition is connected with the rotation of (7r) 
in the spin space, while its absolute value is only slightly changed with temperature. The 
dominant excitations here are the "spin-wave" excitations, except for a narrow critical region 
close to T c . At intermediate values of Q both effects, the rotation of (7r) and temperature 
variation of its absolute value, are important. Thus the l/S'-expansion can be applied to 
describe the temperature phase transition only near QPT. 

Being rewritten through the quantum-renormalized ground-state parameters, finite- 
temperature properties near QPT are universal. The finite-temperature order parameter 
and energy gaps obey the scaling laws 
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S 2 (t_,T) = S 2 (t„,0)f(^-) (55a) 



(55b) 



x cA _ 

A±(t ±J T) = A±(t±,0)g ± ^-£-J 

where A± = A±(t, 0), / and g± are universal scaling functions with /(0) = <7±(0) = 1. The 
transition temperature is determined by zero of the function / (T/cA_o) or, equivalently, of 
<7_ (T/cA_ ) . As discussed in Ref. [p6fl , the functions and <7_(a;) are connected by the 

procedure of analytical continuation. Due to universality of scaling functions (|55a|) , (|55b|) , 
the continuum limit of developed theory, i.e. the action ([34]) , can be used when treating the 
finite-temperature properties. 

Now we pass to calculation of the functions / and g±. Consider first the perturbation 
approach. We obtain 

Al(t + ,T) = A 2 + (i + ,0) + ± f'F, (56) 

for the disordered phase and 

^ 2 , _ N _a , ^ u {2T\ d - 1 „ ( c 2 d 



SUt-,T) = ~ j- d [—) F d -L (57a) 



Alit^T) = s z R (t~,T) 

for the ordered phase, where 



A 2 (t„,0) 3n/2T\^ /c^_ 
^(t_,0) 4lcJ ^21**" 



(57b) 



d— 1 ,7 

= Kd / ( C ° th + * ~ ' F3(0) = 24 

o v y 

and F' d (x) is the derivative with respect to x. Thus we have for the static susceptibility in 
the disordered phase at IA + < T < / 

X " = 7oA 2 ,(t,0)+7(2T/c) d - 1 (58) 

with 7 = 3/o-Pd (0) /4c. 

At T > / thermodynamic properties cannot be determined correctly from the above 
approach since in this temperature region higher-order terms in the 1 / 5-expansion contribute 



to partition function and such an expansion becomes inapplicable. However, one can expect 
that at T ^> / the thermodynamics is the same as for the well-studied Ising model. In 
particular, the susceptibility obeys the Curie law 

X~ = (59) 

on both sides of QPT. 

The equation for the transition temperature reads 

Since at small t_ and d < 3 one has from (^2]) S 2 (t^, 0) oc (f_)^ -1 ^ 2 , we obtain 

T c oc yjtl, 1 < d < 3 (61) 



where the coefficient of proportionality is determined by the solution of Eq. (|60|) . For d > 3 
we derive 

T c (x{t^f'^ l \ d>3 (62) 

The mean-field logarithmic behavior (|5|) is reproduced only for d — *■ oo. 

Consider now the renormalization of the finite-temperature properties at d < 3. To 
this end we use the approach of Ref . |26j , which treats the renormalization of the effective 



classical model. The disordered phase was considered in details in Ref. p6 . Instead of 



analytical continuation of these results to ordered phase, we perform direct calculation of 
finite-temperature properties in the ordered phase. This gives a possibility to calculate cor- 
rectly the value of T c not too close to QPT and also to describe finite-temperature properties 
at T < T c . The generalization of the approach of Ref. |26| to the ordered phase is trivial. 



We integrate out all the modes with nonzero Matsubara frequencies from finite-temperature 
partition function to obtain the effective action for the field 

rc/T 

n= / dT7l x , (63) 

Jo 
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which corresponds to u n = mode, in the form 



S cl = — / d d r 
2T 



K(VU) 2 + Rtf 



3! T 



+ ZiT /,/ ' ,, ' fll + - 



(64) 



Here II = II — II, II is determined by the condition of absence in ( |64"D of terms that are linear 
in II. The parameters of the model (|64]) are given by 

(65) 
and for d = 3 



R(T)= l -U(T)U 2 (T) 



n 2 (T) 



U(T) 



u 



87r 2 g R 



1 8n 2 g R 

3 1 + (3g R /2) ln(l/A_ ) 



T 



2 * ^AV^ 



l + (3W2)ln(l/A_ ) 



1 + 



\ 2T 2 J 
~ /3A 2 r' 



l + (3W2)ln(l/A ) d V 2T 2 



(66) 
(67) 



where 



F 3 (x) = K 3 / q 2 dq 



coth y/q 2 + 



x 



\Jq 2 + 5 



g 2 + x g 2 



(68) 



Fg(x) means the derivative with respect to x, and we have represented (|66D and (|67|) in the 
scaling form by replacing £_ — > A 2 ^ in arguments of F 3 (x), F 3 (x). Near QPT (i.e. for small 
A_ ) function R(T) coincides with that determined by continuation from paramagnetic 
phase, as it should be. The value of K(T) will be needed only in zeroth-loop order, K — 1. 

The critical temperature is determined by the condition H(T C ) = 0. Closely enough to 
the critical point (at ln(l/A_o) 3> 1) we have 

3 



2tt 



(69) 



in agreement with Ref. p6j (our definition of A_ differs (2c?) 1//2 c times from that used in 



Ref. [26]). At the same time, the expansion in the bare splitting (magnetic field) yields 



T c oc c^/tlln 1/3 (l/t_) oc S R (t_,0) 



(70) 



20 



where the coefficient of proportionality can be determined numerically from ([45]) and (|69D. 
Thus, due to ground-state renormalizations, the dependences of T c on the bare and renor- 
malized splittings turn out to be different. 



The resulting classical action (11041) is renormalized in a standard way p3[. One can 



introduce the renormalization constants for finite-temperature theory in the form 

R = (Zj/Z T )R r , IT = Z T Tl r , U = -/-—{Zl/Z T2 )U r (71) 

where the index "r" stands for the quantities renormalized by temperature fluctuations, and 
e = l + e. The expressions for Z-factors are the same as for the ground-state renormalization 
factors (|40|) with the replacement e — > e. Formulas of RG transformation also have the same 
form (|4iD as for the ground-state properties with t — > R, g — > U etc. However, now already 
at d = 3 (e = 0) we have e = 1 and thus the e-expansion can be used only approximately. 
For the energy gap we obtain in this way the expression 



Al(t-,T) 



R(T) 
6 



$TK Z LJJ{T) 
2cR 1 / 2 (T) 



-1/3 



(72) 



where we have put e = 1. For the temperature-dependent magnetization we obtain 

(73) 



GU(T) 



l + 3TK 3 L 1 U(Ty 2/A 



2cR 1 / 2 (T) 



The values of the temperature-transition critical exponents, 

#r = l/3, i^ = 7/12, (74) 

coincide with those of 2D quantum phase transition (|50|). 

The calculated dependence T C (Q) is shown and compared with the mean- field and HTSE 
results in Fig. 3. One can see that near QPT the dependence T C (Q) calculated from (|66| ) is 
in excellent agreement with HTSE data. At the same time, far from QPT our approach gives 
much larger values of T c , as discussed in the beginning of the present section. The inflection 
point of the curve T C (Q), Q* = 0.35/o, may be approximately related to the transverse-field 
value where the "non-spin-wave" excitations becomes important for description of finite- 
temperature properties. 
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For d = 2 the system is far from its upper critical dimensionality (e = 2) and e-expansion 
becomes inapplicable. Therefore we can perform only ground-state renormalizations in the 
results of perturbation theory (|57a ) and ( |57b| ). In this case the critical exponents of the tem- 



perature phase transition still have their Gaussian values. However, universality hypothesis 
predicts that the temperature phase transition critical exponents coincide with those for the 
2D Ising model, 

fa = 1/8, v T = 1. (75) 

With account of the ground-state renormalizations, the result for the critical tempera- 
ture near the quantum phase transition (l/A_o ^> 1) takes the form 

T c oc A_ (76) 

while 

T c oc (t_) 5/12 (77) 

in terms of bare splitting (or external transverse magnetic field). The correct description of 
termodynamics below T c in the 2D case is still an open problem. 



III. THE HEISENBERG MODEL WITH STRONG EASY-PLANE ANISOTROPY 

A. Ground-state properties 

We start from the general Hamiltonian of a spin system in crystal field which induces 
the single-site anisotropy, 

W = Kf-?£ J < J i (78) 

z m 

where V c { is the crystal field potential, J are momentum operators, X is the exchange integral, 
and the direction of spin alignment will be supposed along z-axis. In this Section we consider 
the single-site easy-plane anisotropy which corresponds to 
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V ci = Dj2(Jn 2 (79) 

i 

where D > is the anisotropy parameter. For integer values of J the lowest level is singlet. 
In this case with increasing D the model ( |79|) demonstrates at some value D c a second-order 
phase transition from the phase with collinear ferromagnetic order (J z ) ^ to disordered 
phase. At the same time, the quadrupole order parameter 

Q = 3((J X ) 2 )- J(J+1) (80) 

is nonzero in both the phases. For half-integer values of J, such transition is absent since 
the lowest state is two-fold degenerate. In the classical limit J —>■ oo with J being inte- 
ger, we have D c ~ J (J + 1)1 — > oo, so that integer and half-integer values of J become 
indist inguishable . 

For integer J the ground state is \A) = |0), and first excited state is doublet \Bi) = |1), 
\B 2 ) = | — T) where \M) are the eigenstates of J x . For J = 1 passing to the eigenstates \M) 
of J z yields 

|A) = i=(|l) - | - 1», 

|Si> = |0), W = ^(|l> + |-1» (81) 

To consider the vicinity of QPT, we have to generalize the theory developed in previous 
Section on the singlet-doublet case. Further we restrict ourselves to the case J = 1. In the 



initial spin space, QPT in the model (|78|) with (79) is not of orientational character: the 
spins always lie in the easy-plane. Thus the spin-wave theory in its standard form cannot 
properly describe the model ( |78|) near such a transition (see e.g. Ref. 0]). However, as 



discussed in Refs. pO| , |2T|| , this transition can be viewed as orientational one in the complete 
SU(3) space which includes the SU(2) spin subspace. Most convenient way to consider the 
rotations in the extended SU(3) space is to rewrite the Hamiltonian ( [7S|) with the crystal 



field ([79]) in terms of the Hubbard operators X™ 11 = \mi)(rii 
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~% E [W + Xr){Xf + XJ™) + {X? - Xr^)(X? - X/ 1 -- 1 )] (82) 
z <ii> 

The rotation through "angle" 9 in 577(3) space is performed by the unitary transformation 
operator U(9) (see Appendix C). 

Following to strategy described in Sect. 2, we define the ground-state critical value of D c 
from the condition sin 9 = 1 which yields 



2Xn 



1 3A _6J +J k + J^/Jo 
z k 



2£S 



(83) 



where = 2yXo(Xo — Xk), A(= 1) is the formal expansion parameter. The critical value 
obtained from ( jS3|) in the 3D case is D c /2Xq = 0.73, which coincides with the result of HTSE 
19fl . For the ground-state magnetization we obtain 



{J%=o = t-[l + B(t_)] 



B(t. 



A 



E 



2X + (2 - r/ 2 )X k 



kot 



; 1 + V )I Q - X k + V ll/1 6X + x k + x 2 /J 



Si 



(84) 



(85) 
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where the excitation spectrum is given by (|C12|) , ( |C13|) , and t- = 1 — (D/D 

rj = (1 — i-) 1 / 2 The excitations of a-type have a gap; they are analogous to the excitations 
in the transverse-field Ising model, considered in previous section. The excitations of /3-type 
are gapless due to spontaneous breaking of rotational symmetry in the y-z plane of spin 
space; these excitations are specific for n > 2 systems. Near QPT we have E±p w E^ and we 
return to the perturbation result (|B19| ) for the transverse-field Ising model with E^ = Ek a 
being the critical mode. However, the renormalization of (|85| ) is performed in a different 
way in comparison with the transverse-field Ising model because of another symmetry of the 
model (see below). The energy gap A_ in the ordered phase is given by 



A*(t_,0)=t_ [1+A{t.)] } 
A{t^) = ^-(A + A l )-1 



(86) 
(87) 



where Aq,i are determined by (|C17|) and (|C19| ). 
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In the presence of longitudinal magnetic field, i.e. of the field H, directed along 2-axis, 
both modes a and (3 becomes gapped, since this field breaks the rotational symmetry. As 
well as for the transverse-field Ising model, we can expect that at the intermediate magnetic 
field values 1 — D/D c <C (H/D c ) 2 ^ 3 -C 1 the ground-state properties near QPT will be 
determined by the magnetic field rather than by t±. In this region we obtain the energy 
spectra 



E 



ktv 



4Jn 



i {i + e 2 H )-i k {i-e 2 H /2) 



2(T - J k ) + 0|(2o + J k )/2 2J + e 2 H (Io - J k )/2 



(88) 



where 6h = {4Ji j 'D c ) 1 ^ . Performing the calculations which are similar to those for the 
transverse-field Ising model, we obtain the result 



s = e H 



1 + B\Bl) 



(89) 



where 



3^ k 



2(T +T k )(l + ^/2)-T k (l-^) 



Eh 



ka 



(2 + e 2 H /2)i - j k + (i - e 2 H /2)i 2 /i 



Eh 



61 +l k +l 2 /l 

K 



(90) 



A more complicated situation takes place in the case of the transverse field directed along 
x-axis |^U|,^8J. This field induces a deviation of spins from easy plane. With increasing the 
field value there occurs a cascade of J second-order phase transitions from ferromagnetically 
ordered phases with (J z ) ^ 0, (J x ) ^ to phases which are ordered only along x-axis 
((J z ) = 0, (J x ) 0) and vice versa. The reason for this is the modification of level scheme 
in the magnetic field directed along the hard axis: in the case where lowest state is doublet 
the long-range order along z-axis is present, while in the case of singlet ground state it is 
evidently absent. We do not consider these transitions here (see discussion of such transition 
in Refs. QHIIl). 
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B. Ground-state renormalizations 



The above theory can be easily reformulated in the path integral formalism. The partition 
function has the form 

Z = J D[a,a\b,b^exp^a^ + b^-H(a,a\b,b^ (91) 

where 7i(a,a\b,b^) is the average of the boson Hamiltonian over the coherent states \a,b) 
|JI| (see also Ref. P9fl ). The continuum limit of the theory can be obtained if we introduce 
real variables TT x , y and Q x ,y instead of the complex ones a, b by the relations 

a = tt x + iQ x 

b = Q y + m y (92) 

(Note that n x and ir y correspond to S z and S y in the original spin space, and two additional 
variables Q x>y arise due to passing from SU{2) to 577(3) space). We obtain 

1 C/T 

S C ont = - I dr I d d r 



2iQ x (dn x /dr) + 2iQ y (d7T y /dr) 
o 

c/T 



+Q 2 + (Vtt-) 2 + m 2 7f 2 ] + | J dr J d d rn 4 

o 

c/T 

+h J dr J d d r tt x (93) 
o 

where we have introduced the notations 7r 2 = (l/c )ir 2 , Q 2 = (X /c )Q 2 , h = H / '{Icq) 1 ^ 2 , 
the bare spin- wave velocity is given by Co = 2\p2dT and we have included into the term 
connected with the external magnetic field H along the S z axis. The parameters of this 
model, determined by the continuum limit, read 

m cont = ~~t-d 

Mcont = 6d(c /l)\( 

Ccont = C (94) 
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with £ = 1 — t_ in the ordered phase under consideration. Proceeding in the same way as 
in the previous Section we integrate over Q x>y . Then we obtain the action of the standard 
two- component quantum 4 -theory in an external field, 

c/T 



S cont = l -JdrJd d r [(dn) 2 + m 2 * 2 
o 

c/T c/T 



There is a crucial difference from the one-component model of the previous section, which 
is due to existence in the ordered phase of the gapless Goldstone mode at H = 0. This mode 



changes the renormalization conditions since it leads to infrared divergences [[37]]. To treat 
these divergences, we take the value of magnetic field H finite, but small enough to satisfy 
(H/D c ) 2 / 3 <C 1 — D/D c . The renormalization of the action ( p5|) is considered in Appendix 
D. We obtain for effective Hamiltonian parameters at the scale Ap, A = (2c?) 1 / 2 the results 
{d = 3) 

9? = 9r [1 + (3W4) In(lA-) + (W6) H^/p)) 
V 1 = *k I 1 + (W 4 ) In(lA-) + (to/6) ln(l/p)] 

x [1 + (5W6) ln(lA-)]" 3/5 *o(g R , ti) (96) 

where the function Q (g,x) is given by ( |1J6| ) , 

9r = K^Z^Ucont (97) 

is the renormalized coupling constant. For the non- universal Z-factors we have 

Z L = l 

(ZZ,)- 1 = l + l(i_) + flnl 

(ZI,)- 1 = 1 + A(tJ) - B{tJ) + | ml (98) 

where A(£_), B{tJ) are given by (§^), (|8"T|), the tilde sign means that the contributions of 
the Goldstone mode j3 should be excluded from A{tJ). 
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Putting in the above expressions p = h 1 ^ 2 where h = H/{IqJr(H = 0)] we have for the 
magnetization at d — 3 the result 

3/10 



Mt-,0) 



1 + 

6 t- 



(99) 



(the terms divergent in H are canceled in Jr). The gap for a-type excitations, which 
determines the longitudinal susceptibility, reads 



A 2 (t_,0) = 12Q 2 /\n(l/h) 
t- 



(100a) 



G 2 - 

u °" Z 



L2 



1 + ^lni 

6 t- 



3/5 



9r 



-f2 



J L4~ 



9r 



Up to some nonuniversal factor Z p we have in the one-loop order 0o = Z p (p s /6dcY^ 2 with 
p s being the ground-state spin stiffness. At H — > the gap vanishes as ln _1 (Xo/if) , which 
is a consequence of degeneracy of the system. 

For intermediate values of external mag netic field, i.e. at 1 -D/D c < (H/D c f' z < 1, a 
characteristic scale for both types of excitations is 1/6h, and the expressions for renormal- 



ization factors Zf have the form, which is standard in the two-component (jf model [33 



Then we obtain for the magnetization 



Jr(&h, 0) = 9 H \ 



'^4 



1 + (5W6) ln(l/^) 



3/10 



(101) 



with 



Z L jZ L2 = l-B'(ej I )+g\n— , 



(102) 



7 H 



C. Finite-temperature properties 

Using perturbation theory we obtain for the finite-temperature magnetization the result 
(see Appendix C) 
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{J Z Y = {J Z ) 2 T=0 



To A (IT 



2T\ 

2Cq \ Cq J 



d-i r 



eld 

2T 2 



t + F d (Q) 



(103) 



where £_ = 1 — (D/D c ) 2 . The first and second terms in the square brackets correspond 
to contributions of a and /3-type excitations respectively. At extremely low temperatures 
T <C Zo£- the contribution of a-excitations is exponentially small and the temperature- 
dependent part of magnetization is determined entirely by second term in square brackets 
of ( |103|) . At temperatures T 3> Z t_ the situation changes and both type of excitations give 
the same temperature dependence, the contribution of the mode a being three times larger. 
Consider now the renormalization of the finite-temperature theory. Integrating out the 



field n(q, u n ) with u n ^ from the action (|9~5|) we obtain the action of effective classical 
model 



>cl 



c 

2T 



E g 2 n q n_ q + 



R(T) + 



3h 
W) 



n(T) 



(104) 



ricf/ 



E (n qi n q2 )n 1 , q3 ,5(q 1 + q 2 + q3) + ... 



qiq2Q3 



E (n qi n q2 )(n q3 n q4 )5( qi + q2 + Q3 + q 4 ) + ••• 

qiq2q3q4 

where the field II = 11+ (II, 0) is now two-component one, and the dots stand for higher-order 
terms. For the parameters of the model ( |104j ) we have 



R(T) = V(T)f/(T) 



(105) 



and for d = 3 



n 2 (T) : 189R 



u 



B 2 



U(T) 



8tt 2 ^ 
\n(l/h) 



9 



20vr 2 ^ 



'h 2 c 2 ' 
AT 2 



(106) 



31n(l//i) 

Note that both R(T) and £7(T) vanish at H — > as ln _1 (Xo/i?) due to quantum fluctuations, 
while II(T) is finite in this limit. The value of T c , as determined by the condition II(T C ) = 0, 
read 

3\/3 



2tt 



-c0 r 



(107) 
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where Go is given by fl87P . Thus the result ( |1 7| ) coincides with that obtained in Ref. 
|26|j up to the nonuniversal factor Z p As well as for the 3D transverse-field Ising model, in 
one-loop order the transition temperature turns out to be proportional to the ground-state 
magnetization. 

It should be noted that, owing to presence of the gapless Goldstone mode, the model 
( |104|) is applicable at T = T c only very close to QPT, unlike the corresponding model ( |64|) 
for the one-component case. Thus one can put Z^ = Zn = 1- The calculated dependence 
T C (D) is shown and compared with HTSE data in Fig. 4. As well as for the one-component 
case, the result (|107|) agrees well with HTSE data closely enough to QPT. A more complete 
treatment can be performed by considering quasimomentum-dependent vertices in ( |104|) . 
This is a complicated task which is not considered in the present paper. 

By the same reason, the model ( |104| ) cannot be used for determining magnetization 
below T c . However, one can expect from the hypothesis of universality the standard value 
of the two-component three-dimensional 4 theory critical exponent 

Pt = ; 6 s = 7/20, (108) 

H 2 2(n + 8) 1 ' V ' 

which is practically the same as in the one-component case (/?t = 1/3)- Unlike the one- 
component case, the logarithmic correction to S(t_,T) is expected near the temperature 
phase transition due to the gapless Goldstone mode. 

For d = 2 the contribution of the gapless mode is logarithmically divergent and therefore 
the long-range order at finite temperatures is absent, unlike the case of the transverse-field 
Ising model. 

IV. CONCLUSIONS 

In the present paper we have considered systems that demonstrate in the ground state 
a quantum phase transition (QPT). Near QPT the saturation moment Sq is small, but the 
Curie constant in ( [5§D is not suppressed. We have T c oc Sq which is determined by the 
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value of dynamical critical exponent, z = 1. The susceptibility ( |58| ) in the intermediate 
temperature region A <C T/I <C 1 is determined by the small ground-state energy gap 
Ao and demonstrates a 1/T d_1 behavior. In the strong enough longitudinal magnetic field 
Ao <C (H/IS) 1 / 3 <C 1 the ground state parameters are determined by magnetic field value 
rather than by closeness to QPT. The corresponding dependences have been obtained. Our 
approach gives a possibility to investigate both non-universal and universal renormalizations 
of the ground-state parameters. The ground-state renormalizations turn out to be important 
in the vicinity of QPT. Thus the results for thermodynamic quantities (e.g. transition 
temperature) have different forms as functions of renormalized splitting and bare transverse 
(external magnetic) field, see (|69D,(|70D and (17(f), (fTTj). This should be taken into account 
when treating experimental data. 

The discussed class of magnets is similar in some respects to weak itinerant magnets. 
Note that for weak itinerant ferromagnets we have T c oc S (see, e.g., Ref. p2[), which 
is due to that main contribution to thermodynamics comes from paramagnons (z = 3). 
As well as for considered localized-moment systems, calculation of non-universal ground- 
state parameters for itinerant magnets is of interest, in particular for different forms of bare 
electron density of states. 

Now we discuss the experimental situation for some systems exhibiting magnetic and 
structural transitions. The compound DyVC>4 demonstrates a structural phase transition 
at Tjj = 14 K. The low-lying energy levels in the spectrum of this system are two Kramers 
doublets with the splitting Ai = 27cm~ 1 at T = and A 2 = 9cm _1 at T > Tp. Neglecting 
the Kramers degeneracy one can describe this system by the transverse-field Ising model 
with Q/I = 1/3 (see Ref. JF]). The corresponding point in Q/I-T c coordinates is marked in 
Fig. Id. This point lies exactly on the HTSE curve and therefore HTSE results are applicable 
in this region of Q/I. One can see that DyV04 lies far from QPT, so that the above-developed 
theory is not applicable for this system. 

Other systems, which are well described by the transverse-field Ising model, are the fer- 
roelectric quantum crystals like KH2PO4 (see, e.g., Ref. f|). However, to our knowledge, 
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corresponding detailed data on ground-state order parameters are absent. To fit experimen- 
tal data on T c of Ref. 0], we need explicit dependence of the tunneling parameter Q on 
pressure. 

There are very few experimental data on singlet-singlet systems demonstrating magnetic 
phase transitions. The system LiTb :r Y 1 _ :r F4 EDI is usually assumed to be characterized 



by long-range exchange interactions and therefore well described by the mean-field theory. 
The singlet-doublet case is represented by the system NiSi2F 6 which is an J = 1 easy- 
plane Heisenberg magnet. The anisotropy constant is changed under pressure and thus the 
value of D/D c can be varied near unity in the experiment. The pressure dependence of 
anisotropy constant was measured experimentally @. However, to our knowledge, the data 
on the pressure dependence of exchange parameter are absent, although it is supposed to be 
considerable ]n|. There are also few experimental data on the ground-state magnetization 
near QPT. At p — 8.6Kbar, one has the experimental values J(T = 0) = 0.3 and T c = 
HOmK |2j. The calculation according to ( |107| ) yields I = 60mK which is larger than the 
p = value, X = 40mK |42|| . Praseodymium in the dhcp phase contains both "cubic" and 
"hexagonal" sites flip^l, so that separation of different contributions makes an additional 
problem. Generalization of our approach to the singlet-triplet case in connection with the 
Pr ions in cubic crystal field will be presented elsewhere. 

Generally, the 1/5* perturbation theory combined with field-theoretical scaling analysis 
enables one to obtain a description of ground-state properties of transverse-field Ising model, 
which in a good agreement with the results of the fourth-order ground-state perturbation 
theory |13| for all the values of Q. The only fitting parameter used is the critical field value fl c . 
The finite-temperature properties are considered with the use of approach of Ref. |26) . The 



same analysis for the 5 — 1 easy-plane Heisenberg model is performed within the expansion 
in formal parameter A(= 1) which plays the role of 1/5. In this case, besides the critical 
mode, there is a gapless Goldstone mode, which considerably modifies the conditions of 
renormalizations. The consideration of QPT in degenerate systems induced by the external 
magnetic field within the approach used is of interest. In particular, in the case of single-site 
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anisotropy oscillations of the effective moment with increasing magnetic field or temperature 
are expected in such systems with S > 1. 

It is of interest to apply the approach used to various 3D and 2D systems demonstrating 
orientational and metamagnetic phase transition with changing the external magnetic field 
or anisotropy, e.g. for yttrium garnets [|J and magnetic films [[0]. Depending on a concrete 



physical situation, such systems can be described by the strongly anisotropic transverse-field 
Ising model or Heisenberg model with small anisotropy. 

APPENDIX A: MAPPING OF THE ANISOTROPIC HEISENBERG MODEL 
ONTO THE TRANSVERSE-FIELD ISING MODEL 

In this Appendix we discuss a possibility of mapping procedure of anisotropic Heisenberg 
model ( |78| ) onto the transverse-field Ising model. We consider only one important case where 
the lowest level of V c f is singlet and there is QPT to disordered phase at strong enough V c f. 
Provided that the first excited state is also singlet, neglecting all energy levels except lowest 
and first excited states we can introduce the pseudospin-1/2 operators S, to obtain PIJlQfl 



Kf=-A]T^, J* = 2aS* (Al) 

i 

where a = (A\J Z \B) is the matrix element of J, \A) and \B) are the lowest and first excited 
states, A is the energy gap between these states. (It should be noted that left-hand sides 
of Eqs. (|A1|) act in real-spin space, while right-hand sides in pseudospin space. Thus 
the equality signs are used only in the sense of identity of averages.) Then we obtain the 
transverse-field Ising model with / = 4a 2 X and Q = A. The order parameter of Heisenberg 
model (J z ) is connected with the order parameter in the transverse-field Ising model by 

(J z ) = 2a(S x ) (A2) 

Consider now the case where the excited state is a multiplet with the states \B m ), m = 
1...N — 1. Neglecting the degeneracy of upper energy level, one can use the same mapping 
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( |Al| ) if we choose a 2 = Ylm=i(A\J z \B m ) 2 . However, in this case the original SU(N) spin 
space is projected onto SU(2) pseudospin space, and thus N — 2 degrees of freedom are 
neglected. Thus this approach does not give a possibility to take into account properly the 
symmetry of the original model and therefore can be applied only outside the critical region. 
To obtain a correct description of such systems in the critical region one should consider the 
transition in complete SU(N) space. 

The above consideration shows that, in principle, the transverse-field Ising model ([I]) can 
qualitatively describe singlet magnets even in the case where the exchange interactions in 
the true momentum space are isotropic, as in model (|78|) . 

APPENDIX B: CALCULATION OF SPIN GREEN'S FUNCTION AND ORDER 
PARAMETER OF THE TRANSVERSE-FIELD ISING MODEL WITHIN THE 1/5 

EXPANSION 



r 2 . 

XI 



Consider first the disordered phase where (ir x ) = 0. Representing 7T zi = (1 + 1/S — tx\ 
Kyi) 1 ! 2 and assuming {n 2 ^) ~ 1/S (validity of this statement will be checked below) we 
expand square root to second order in 1/5* to obtain 



l/T _ l/T 



i$ w Sy* f 9-Kyi _ dir xi \ t iS ^ [ , ji \ ( „ d^vi <- )77 ->-> 

* " ' ' I. 



s dyn = E Mm ~ ^i^r + t £ / dr vli + 7r - 



dr nyi dr )' 8 4- J "" y >xi ' dr ^ dr 

o 

l/T 



Set = ~\ i (,T 







IS 2 £ * n * xj + (TV - QSw s ) £ {<, + 4) - X ? + <) 2 



where % = (1 + 1/2S 1 ) -1 and P = XL„ 1 ( w n being the Matsubara frequencies) is formally 
divergent quantity which comes from the measure of integration, this divergence will be 
canceled in final results To first order in 1/S (we suppose that Q c ~ loS") we obtain by 
standard perturbation theory methods the matrix two-point vertex function of n x , ir y fields 
in the form 

' ns(w s + ^)+sr-Tv-i< l s 2 iSuJn (w s + ^f) ^ 

K iSu n (w s + £±*) nS (w s + ^) + ST - TV ) 
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T(q, UJ n ) 



(B3) 



where 



and 



T = i(7r xi (dn yi /dr)) = I £ ~TTm> (B5) 

^ q,^ n + "^q 



with = y — S7q) being the bare "spin-wave" spectrum in the disordered phase. (One 
can easily verify that X, Y are of the order of 1/S, as it was supposed in the beginning) . 
Using the identity 

QSX + ST - TV = I S 2 X' (B6) 

where 

* =<*-■< >=^£^ (B7) 

to eliminate the divergences, we derive 

/ OcL, i X+y\ , ro2y/_r c2 „• o, , /'„„_ , X+Y\ \ 



T(q, w„) 



B8) 



QS {w s + + I S 2 X' - I^S 2 iSu n {w s + 
[ isu n {w s + ns {w S + + I S 2 X' 

The value of Q c is determined by the condition 1^(0, 0) = However, the averages ( p4|) and 



/ 



( [B7D which determine l/S'-corrections to r(q, uj n ) are fi-dependent itself. For consistency, 



we have to calculate these averages with the zeroth-order value Q c = IqS. Then we obtain 
the result for Q c (|l"5|) of main text. At an arbitrary Q we perform the replacement Q — ► 
IqS(Q/Q c ) in ( |B4| ) and (|B7|), which changes r(q, o> n ) in the order of 1/S* 2 only and gives a 



possibility to take into account consistently the shift of Q c owing to quantum fluctuations. 
Substituting the result for Q c into (p8|), we obtain the results (0), ( |16D and fllTD of the 
main text. 
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In the ordered phase we rotate the coordinate system through the angle 9 around 7r y -axis 
and again, in new coordinates, expand tz z in powers of tt x , n y (ify = n y ). Then we obtain to 
fourth order 



IS 2 



l/T 



E / ^ r -Ti^Wg 1 sin 29 + 7i xi 7i xj cos 2 9 - (ii 2 xi + 7r^-) sin 2 9 



<*?> o 



i ~ i~2 ~2 \ • o/i i \ xi ' yiJ\ xj ' yjJ V xi 1 yil ■ 2 n 

+ 2 7T *i( 7l xi + n yi) sm 29 + sm 



l/T 



-05^ / rfr 



7r T ,- sin 



4 



7T 2 - + 7T 2 - (7T Z - + 7T , 

< "yi n \"xi ' yU ri 

-Mis COS 9 — - COS 9 



l/T 

^— E J dT + 



7T 



(B9) 



5d yn having the same form ( pT|) as in the disordered phase with the replacement n — ► 7r. 
Determining the angle # from the condition (vf^.) = we obtain 



cos 9 



ft 



1 + 



1 X + 2X' + Y 
2S 2 



-i 



(BIO) 



where 



* = = T E 



-X"' = (ftxiftxj) = t E 



q,a; n n q 



^ = «)=^E 



q,u; n 



(Bll) 



-E q = S^J Iq(Iq — rjlq) and 77 = (f2// 5') 2 - Besides the averages ( Bll|) , we introduce the 
quantity 



T = i(n xi {dn yi /dT)) = - 



For the two-point vertex function of the fields n x , 7r y we have 
/ 

r(q,w„) = 



/o^ 2 (1 + X - r]X') - I„S 2 W + T qn + ST - TV 



V 



iSu n (w s + 



Y+X 



(B12) 



lSu n (w S + ^) 

I S 2 (1 + Y - r)X') + ST - TP 

(B13) 
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where 



W = 7]{w s + X + 2X' + Y) + (1- rj)X', 



(B14) 



The term with 



S 



F m = -^(1 - V )J2 [(2/k/k+q + 4/ q / k + 2/2 +q + I 2 q )M xxxx (k, q) 



k,ui m 



+/^M rara (fc, q) + 2/ q (/ q + 2I k )M xyxy (k, q) 



(B15) 



where = (k,u m ), q = (q,,U) n ) and 



Mapys(k,q) = (nMn^-k^in^k + q)n s (-k - q)) (B16) 

(a, P,j,5 = x, y) arises due to the contribution of the cubic term in (p9| ) in the second order 
of perturbation theory. This term has the same order in 1/5* as other terms and should be 
retained. Using the identity 



I S 2 X + ST - TV = I S 2 7]X' 



(B17) 



which is an analog of ( |B6l) for the ordered phase, we obtain 
r(q,Wr») 



( S 2 (I -WQ+F m iSu n (w s + X/2 + Y/2) 



V 



(B18) 



/ 



iSuj n (w s + X/2 + Y/2) I S 2 (l + Y-X) 

Performing again the replacement Q —>■ IqS(Q/Q c ) in ([Bll ) and reexpressing ( |B18| ) in terms 
of Q c , we obtain the results (pTT]) , (|19|) and of the main text. For the temperature- 
dependent order parameter we obtain 

S = S(ir x ) = S sin6>(7f x ) 



77 ^ 2/ + / q 



1/2 



+— y 



+ 25 2 



(B19) 



Note that near QPT the last multiplier in this expression is close to unity and only slightly 
temperature-dependent. Thus it can be replaced by its zero-temperature value. 
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APPENDIX C: ROTATION IN SU(3) SPACE FOR THE EASY-PLANE 

HEISENBERG MODEL 



Following to Refs. 
mation 



2TT we perform in the Hamiltonian (RI3) the unitary transfor- 



Xf q = U\6)Xf q U{8) 



(CI) 



with 



U(0)=exp[e(X- 1 > 1 -X 1 >- 1 )/2] 

= 1 + [cos(0/2) - l^X- 1 '- 1 + X 1 - 1 ) + sin(^/2)(X- 1 - 1 - X 1 -" 1 ) 

Then the Hamiltonian takes the form 



(C2) 



(C3) 



with 



1-L {1) = - E [( D sin 9 + 2 J cos 2 9) (2X, 11 + X°°) + DXf° + (D - 2 J sin 6) cos + Xr 1 ' 1 ) 



= -f E [cos0(2X 4 n + X°°) - sin^X/'" 1 + X^ 1 ' 1 )] [cos6{2Xf + X°°) - sin^X, 1 '- 1 + X7 1 ' 1 



<«j> 



J 



-EE [2X° ff XJ° + a sin 9{Xf'Xf I + XfXf) + cos ^X^X" + X '^; 



<ij> cr=±l 

where Xo = 2dX and we have dropped the tilde sign at the X-operators. Further we represent 
Hubbard operators via boson ones pO|,pT| 



X 



1,-1 



i\l-Xa ] a-\b ] b) 1 ' 2 



X -- 1 = 6 t (l-Aa t a-A6 t 6) 



1/2 



X 



10 



a\ X 00 = tfb, X 11 = a f a 



(C5) 



where A(= 1) is the parameter introduced to construct perturbation theory (cf. the Holstein- 
Primakoff expansion in the case of a Heisenberg magnet). Then we obtain the Hamiltonian 
of the bosons 
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Tt = Tii + Ti.2 + Ti.3 + 7^4 + ... 



(C6) 



where 



Hx = - cos 9(D - 2J sin 9) £ (a\ + a 4 ) 

ft 2 = (D sin 9 + 2J cos 2 9) £ a}a< + [D(l + sin 0)/2 + X cos 2 0] £ &}&< 

j i 

-- [ sin2 Q ( a i + a i) ( a ] + a j) + 2b l b j ~ sin e ( b Wj + 
2 (ii) 

7i 3 = — cos 6>(/J — 2X sin 6 1 ) E [aj (a|oj + fttft^ + {a\di + 



+Zcos#E [sin^ (2a\ai + b\b^j faj + a^ — [biO^bj + bjbjaj 



(ij) 



H 4 = -- [ sin2 6 ( a l a i + 2b l b i) ( a i a i + 2b ) b i) + 2a \ a o b ) b i 



(ij) 



— cos 2 9 (b\ + fejj (ajaj&j + + h.c.J + cos 29 fajOj&J&j + h.c. 

— (a\ajajaj + aj&j& 3 -a 3 - + h.c.) + COS26 1 (^a|a|aja] + af&j&jaj + h.c 



(C7) 



(C8) 



(C9) 



(CIO) 



and we have omitted terms containing more than four Bose operators. Diagonalizing the 
quadratic part H2 of the Hamiltonian flCTf) we obtain 



(Cll) 



where the spectra of excitations are given by 



E ka = 2y/I (I -r? 2 / k ), 

= y/(l + ri)(I - I k )[I - /k + r/(/ + 4)], 

?7 = sin# = D/D c . The angle # is determined by the condition (e^) 
the formal parameter A we have 

sm9=^-[l-X(R 1 + R 2 )}- 1 



(C12) 
(C13) 

0. To first order in 



(C14) 



kcv 



k ^k/3 ^ 



(C15) 
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The gap in the excitation spectrum for the mode a to first order in A is given by 



where 



4> = £ 



^k=o, a = 2 v /(l + J Bo + J B 1 )J A_ 
A 2 =(l+A + A 1 )(l-r 1 2 ) 



?/ 2 7k - 4 7k -_4 7k(l + V + rf)/{l + v) ~ (V + 1) 



(C16) 



ka 



Eh 



k/3 



8(l + 7 k)-20r ? 2 7 k + llr / 4 7k | 2(1 + 77) - 2 7k + 2/7 3 7 2 

-^ka -Ek/3 

(5/2) (2 - ry 2 7k ) 2 + 4(2 - r^k^Tk + (5/2)(/7 2 7k ) 2 



(C17) 
(C18) 



+ 1 ^2 (^ 2 + l)(l + ^-7k) 2 + ( ^ 2 + l)(r?7k) 2 + 477(1 + 77 - 7 k)(?77k) 
2 k 



ft 3 



1 



5i = ^(i-^ 2 )E 

z k 



4r/ 2 (2 - 77 2 7k ) (1 + 77 - 7k ) 2 - (r7 7k ) s 



E 3 



ft 3 



(C19) 
(C20) 



The mode (3 has Goldstone type and is gapless to arbitrary order in A. 

APPENDIX D: RENORMALIZATION OF THE TWO-COMPONENT ^ MODEL 
WITH SPONTANEOUSLY BROKEN SYMMETRY 



In this Appendix we consider the renormalization of the action (|95 ) in the ordered phase, 
m 2 < 0. Introducing the quantity k 2 = —2m 2 > and performing the shift tt x — > tt x + tt we 
obtain 

\2 1 /,2 1 qM„2 



S = - I (It I d"r 





(dir) 2 + (ac + 3h)n 2 + hir 2 

c/T 



+| J dr J d d r (4tt x + n 2 ) tt 2 



(Dl) 



where h = h/w , W = (3k 2 /u) 1 ^ 2 and 7T = 7Tq(1 + h/ k 2 ) is determined by the requirement 
of absence in the action of terms that are linear in tt x . In these notations the condition of 
smallness of magnetic field is h 1 / 2 <C k (the condition of closeness to QPT k C A is also 
assumed). Under these conditions, the action ( pi|) has two characteristic lengths, 1/h 1 / 2 
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and 1/k. This situation is the same as in the theory of crossover phenomena [ 36 1 . Further we 
follow to Ref. |36| to include exactly the smaller characteristic length into Z-factors. Then 
we obtain to one-loop order 



7Cont 



/7Cont 
Z 2 



^cont 



0(g 2 
9 



1 



, g_ _ ^9_^_ 

2e (1 + K 2 /p 2 y/ 2 6e 3eA £ ' 

3g 1 ^_ 5g_ff_ 

2e (1 + K 2 //i 2 ) £/2 3eA £ ' 



(D2) 



(note that the Ward identities guarantee that the structure of the interaction term is pre- 
served by renormalizations, and one renormalization constant is sufficient to renormalize 



all four-particle vertex functions, see, e.g., Ref. f33|). The flow functions for the effective- 
Hamiltonian parameters are 



/3(g,K/fi) 

j(g,K/ii) 



-eg 



g 



zg 2 



+ 



r 



2 l + K 2 /p 2 6 

i g 



(D3) 



2 1 + n 2 /p 2 6 

Putting in these expressions e = 0, performing the integration and supposing that scaling 
starts at /i > k we obtain the effective-Hamiltonian parameters at the scale pp 



1 1 

g P g 



ln(p 2 + n 2 /p 2 



6 



In p 



K 2 = k 2 exp 



1 

X- 



y 2 + K 2 /p? + 3p', 
l 



(D4) 



1 - (g/Q) hip' - (3g/A) ln(p' 2 + k 2 / p 2 ) _ 

Our plan now is to use these scaling formulas to reach the scale pp ~ h l l 2 <C k. For these 
values of p the formulas (|D4|) are simplified: 



g- 1 = g- 1 [1 - (3g/2) ln( K //x) - (g/Q) hip] 



k- 2 [l-{3g/2)\n{K/p)-{g/6)\np] 



x[l-(5g/3)HK/p)}- 3/5 %(g,K 2 /p 2 ) 



(D5) 
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where $o(g,x) is given by 

i 



1 - :j1 



Ay + x 



y + x 1 - (g/12) \ny- (3g/A) ln(z + y) 

6{x-y) A9(y - x) 



(D6) 



1 - (g/12) Iny- (3g/A) lnx 1 - (5g/6) \ny _ 

and 9(x) is the step function. At extremely small x we have &o(g,x) ~ exp(l/ln 2 x) ~ 1. 
The result ([05]) with $ (g, k 2 / p 2 ) = 1 can be obtained more directly if we perform the 
scaling procedure in two steps: at the first step p ^> k/{i and the flow functions are the 
same as for the two-component isotropic 4 model, while at the second step p <^ n/p and 
the flow functions include only contributions of the Goldstone modes. The scaling formulas 
are joined at p = n/p. However, this procedure does not give a possibility to describe 
correctly the contribution of the crossover region p ~ n/p. Putting in the above results 
p = A = (2c?) 1 / 2 , we obtain the result fl96|) of the main text. 
At finite but small e we obtain in a similar way 



a; 1 = g- 1 



.-2 



K 



-2 



1 + (3g/2e)(K^/p- £ - 1) + (g/6e)(p- £ - 1 
1 + (3g/2e)( K - £ /p- £ - 1) + (g/6e)(p- e - 1 

3/5 



X 



l-(5g/3e)(K- £ /p 



(D7) 



(D8) 



with some function § £ (g,x) ~ 1. 
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One can obtain the same results in a more correct way (which is free from divergences) 
by keeping along the calculation the imaginary time to be discrete, r = r, = i/MT with 
i = 0..JV, and passing to the limit M — > oo only at the end. 

Figure captions 

Fig.l. Ground-state energy gap E (Q) (Q > fl c , a) and order parameter 5(0) (Q < Q c , 
b) for the 3D transverse-field Ising model in different approaches. The value Q c = 2.58/ is 
used for calculating the 1/5" and RG' curves. Arrow shows the value of the critical field Q c , 
obtained by 1/ S'-expansion 

Fig. 2. Ground-state energy gap Eq(Q) (Q > f2 c ,a) and order parameter S(Q) (ft < Q c , 
b) for the 2D transverse-field Ising model in different approaches. The notations used are 
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the same as in Fig.l 

Fig. 3. Transition temperature as a function of Q/I for the 3D transverse-field Ising 
model in different approaches. 

Fig. 4. Transition temperature as a function of D/X for an easy-plane ferromagnet in 
different approaches 
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